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A three-dimensional scalar field theory model of center vortices and its relation to 

k-string tensions 

John M. Cornwall* 
Department of Physics and Astronomy, University of California, Los Angeles CA 90095 

In d — 3 SU(N) gauge theory, we study a scalar field theory model of center vortices, and their 
monopole-like companions called nexuses, that furnishes an approach to the determination of so- 
called fc-string tensions. This model is constructed from string-like quantum solitons introduced 
previously, and exploits the well-known relation between string partition functions and scalar field 
theories in d = 3. A basic feature of the model is that center vortices corresponding to magnetic 
flux J (in units of 2n/N) are composites of J elementary J = 1 constituent vortices that come in 
N — 1 types, with repulsion between like constituents and attraction between unlike constituents. 
The scalar field theory is of a somewhat unusual type, involving N scalar fields <j>i (one of which is 
eliminated) that can merge, dissociate, and recombine while conserving flux mod N . The properties 
of these fields are deduced directly from the corresponding gauge-theory quantum solitons. Every 
vacuum Feynman graph of the theory corresponds to a real-space configuration of center vortices. 
We use qualitative features of this theory based on the vortex action to study the problem of fc-string 
tensions (explicitly at large N, although large N is in no way a restriction on the model in general), 
,—1 , whose solution is far from obvious in center-vortex language. We construct a simplified dynamical 

picture of constituent-vortex merging, dissociation, and recombination, which allows in principle for 
the determination of vortex areal densities and fc-string tensions. This picture involves point-like 
"molecules" made of constituent "atoms" in d — 2 which combine and disassociate dynamically. 
These molecules and atoms are cross-sections of vortices piercing a test plane; the vortices evolve in 
a Euclidean "time" which is the location of the test plane along an axis perpendicular to the plane. 
A simple approximation to the molecular dynamics is compatible with fc-string tensions that are 
^li- . linear in k for k <C N, as naively expected. 

oo ; 

Q PACS numbers: 11.15.-q, 12.38.-t, 11.15.Tk UCLA/04/TEP/5 

o 

^- ' I. INTRODUCTION 

o 

A. General remarks 

Mh, Center vortices have been intensively studied since their introduction more than twenty- five years ago [U 0- LJ- [H S Lj| • 

There are now many lattice studies (for a modern review, see Q), which confirm the basic topological mechanism for 
confinement in center-vortex theory. 

In this paper we present a model for center vortices and nexuses in d = 3, based on the well-known connection [8j 
between closed-string partition functions and (scalar) field theories in three dimensions. This model is applicable, in 
j_j i principle, for all SU(N) gauge theories. In d = 4 there is an analogous model, not explored here, which would be a 
closed-string theory of unusual type. In both d = 3 and d = 4 we would find exactly the same picture described below 
of "atoms" and "molecules" interacting in a plane, as a description of center-vortex areal densities. These densities 
are essential to our understanding of confinement. 

While it is desirable in general to have a good model for center vortices that is simpler than directly solving the full 
gauge theory, the real point of the model is the extent to which it can clarify the full theory. We have in mind here 
applying the model to the question of string tensions with non-trivial JV-ality, which has presented some difficulties 
in center-vortex theory. This remains one of the outstanding first-principle problems of confining gauge theories in 
the center vortex picture. 

The primary focus of the present paper is the model itself, but before presenting the model we will make some 
remarks on the problem of the so-called fc-string tension, which is the string tension a(k) that arises between a test 
particle (in SU(N) gauge theory) whose group representation is described by a Young tableau of a single column of 
k < N boxes, and its antiparticle, described by N — k boxes. This test particle and its antiparticle form a singlet and 
are confined. The fc-string problem exists only for N > 4, since for N = 2, 3 the k representations consist only of the 
fundamental and antifundamcntal representations. But for larger N one expects to have different string tensions for 
different k. 



Email: Cornwall@physics.ucla.edu 



Various simple arguments, for example, standard large- N considerations, suggest (at least if fc <SC N) that a(k) ~ fc. 
Or one could say that a k representation is the totally antisymmetric product of k fundamental representations 
(quarks), so that in terms of the mesonic string tension a (I) the fc-string tension ought to be composed of k quark- 
antiquark strings and hence k times as big. 

Unfortunately, it has not proven to be so easy to derive u(k) ~ k in the center vortex picture of confinement. The 
reason is that in the center vortex picture confinement arises from the topological linking of center vortices with Wilson 
loops, and it is not easy to envisage the k separate quark-antiquark strings mentioned above. In the center vortex 
picture, a quantity from which cr(fc) can be found is the discrete Fourier transform of the areal (two-dimensional) 
density of center vortices with different values L of magnetic flux, measured in units of 2ir/N, and one needs to find 
these areal densities as a function of L. 

With no compelling picture of vortex areal densities, the result so far has been that theoretical approaches to cr(fc) 
in the center vortex picture have been largely phenomenological. Workers often take Casimir scaling of <r(k) or analogs 
to M-theory as the two possible standard forms for the fc-string tensions. These forms are 

k(N-k) smjirk/N) 

vcas(k) = o-i— — , <7 M th(k) = cri . . (1) 

N — 1 sin(7r/iV) 

Both of these lead to er(fc) ~ k for fc <C N, but they have different terms at non-leading order for large N. There are 
arguments 9] that non-leading terms exclude Casimir scaling, but in fact which of these two forms, or some other 
form, actually holds is not known now. In any case, we will not study non-leading terms here. Both forms above 
also obey the necessary group-theoretic conjugation symmetry a(k) — a(N — fc). Certain works that address c(fc) in 
the center vortex picture depend on some form of quadratic approximation that inevitably leads to Casimir scaling 
[l(j ; others are phenomenological, and ask what density of vortices is needed to produce c r(fc) as given by the above 
choices [H|,[l2j or other related forms. Using the dilute gas approximation (DGA), Ref. [IlJ finds unphysical areal 
densities amounting to the unchecked pileup of center vortices at a single point. Ref. [IJ (hereafter GO) supplies a 
relation between fc-string tension and areal densities which is usable in principle where the DGA is not, and we will 
make use of this relation in the present paper. 

In addition to these theoretical works, there are lattice simulations |13l HI . Ha . flq for N = 4, 5, 6. Most of this 
lattice work supports the M-theory form of the fc-string tension, but it does not seem possible to draw unassailable 
conclusions from the present data. A very recent work 17] concludes that lattice data lie somewhere between M-theory 
and Casimir scaling. 

So far (to the author's knowledge) there has been no discussion from the fundamental principles of center-vortex 
theory of the areal densities needed to find the fc-string tension. While we do not claim to solve the fc-string tension 
problem here, we do offer a path toward solution in the present paper. 

B. Outline of the paper 

The main purpose of the paper is to present the d = 3 scalar field theory that models the gauge theory. [There 
are other models of center vortices; see [l8| for SU(2) and [n| for SU (3). But there is no fc-string problem for these 
gauge theories.] A secondary purpose is to introduce the possibility of exploiting the scalar field theory to discover 
something about fc-string tensions, although we do not go very far in this direction, and restrict our considerations 
to large N. [It is not yet known whether center vortices become unimportant at large N, although it is very possible 
|20j that they survive in this limit.] After a quick review of the gauge-theory form of center vortices in sec. [HI we 
show how to express these gauge-theoretic results via a scalar field theory in sec. IIIII This field theory is constructed 
directly from the gauge-theory picture of the quantum solitons constituting center vortices and nexuses 20]. The 
basic principles of equating a theory of closed strings, such as center vortices, to a scalar field theory have long been 
known [8j, but the theory we find is considerably more complex than envisaged in these pioneering references. The 
scalar theory transcribes the properties of vortices (strings in d = 3) and nexuses (points in d — 3) into Euclidean 
world lines for the particles of the scalar field theory. These strings, with nexuses and antinexuses sitting on them, 
undergo processes whose essential features we will incorporate into the scalar field theory. Specifically, every vacuum 
graph of this scalar field theory is in one-to-one correspondence with a process of physical vortex strings merging, 
dissociating, and recombining. Every center vortex of flux J > 1 is describable as a composite of J constituents, each 
of unit flux. 

Sec. IIVI goes into some kinematical simplifications which are important for our specific implementation of the model 
in determining fc-string tensions, but it is better to outline first some more general features of the relation of the scalar- 
field theory model to string tensions, as put forth in sec.0 These general relations will survive the simplifications 
of sec. IIVI To describe string tensions, which involve the piercing of d = 2 surfaces (Wilson loops) by center vortices, 
we introduce certain probabilities p(L), where p(L) is the probability that a center vortex of flux L intersects a large 



but finite planar surface, which we call the test plane. These probabilities are converted to areal densities by dividing 
by a squared correlation length A 2 , as GO describe. We construct a continuous series of cross-sections of d = 3 
vortex processes by rigidly translating the test plane along its normal (say, the z axis). The cross-section of a vortex 
in the plane is point-like, and we refer to composite vortices as molecules and the unit-flux constituents as atoms. 
Because of the underlying random-walk nature of the vortices, these atoms and molecules appear to diffuse in the 
test plane as z changes, from time to time colliding and interacting with each other. Under the assumption that 
such processes are ergodic, averaging over z should yield the same results for areal densities as averaging over many 
gauge configurations. We construct a a simple master rate equation in z, in an attempt to describe these merging, 
recombination, and dissociation interactions. Which processes are allowed and which are forbidden is governed by a 
set of force and action rules derived directly from the actions of the vortices in the process. The master rate equation 
yields (in principle, at least) the vortex areal densities as densities of atoms and molecules in the test plane. 

At this point we do not have enough quantitative information about the scalar field theory model (couplings and 
masses) to justify a thoroughgoing study of the master rate equation. Instead, we will make a few suggestive remarks, 
based on some kinematic simplifications which are motivated by large-TV considerations. These simplifications lead 
to a modified form of the DGA, which is free from the objectionable features of the phenomcnological probabilities 
p(L) found in 11] using the standard DGA approach. 



The kinematic issues are set out in sec. IIVBI following GO. These authors formulate the fc-string question in a 
way that is not subject to the low-density restrictions of the DGA. They go on to give purely phenomcnological areal 
densities which both fit a linear rise of a(k) and are well-behaved in the large-TV limit. They define p(L) as the 
probability that a square of side A is pierced by a vortex of flux L, where A is the correlation length of vortices. [The 
probability p(0) is the probability that the square has no vortex in it at all.] They then argue that 

exp[-a(fc)A 2 ] = f>(i)exp[^]. (2) 

o 

The inverse discrete Fourier transform yields p{L) in terms of o~{k) as 

1 N ~ 1 -9 l T 
PM = ± E ex P[^?F] ex P [- CT (fc)A 2 ]. (3) 



N ^ L N 

fc=0 

Conjugation symmetry applied to this inverse transform then yields p(L) = p(N — L) for the areal densities, and then 
Eq. J2J yields a(k) = a(—k). In fact, as is well-known, a vortex of flux TV — L is actually an anti- vortex of flux —L, so 
we should also have p(L) — p{—V). All these conjugation-symmetry relations will hold in our basic model for center 
vortices. 

We are interested in large N, and in sec. UVBl we explore the consequences of replacing discrete variables and sums 
over these variables by continuous variables and integrals over these variables. In so doing, we lose track to some 
extent of conjugation symmetry, but recover enough of it for our purposes in the evenness of p(L) and a(k). First, 
note that Eq. © suggests a scaling symmetry, in which p{L) is replaced by a continuous function: 

p(L,N)->±p(x), x=^ (4) 

whereby the sum in the fundamental GO equation J2J) is replaced by an integral from to 1 over x. Conjugation 
symmetry for p{x) is retained in the form 

p(x) = p(l-x) = p(-x). (5) 



However, at this point we have lost (see sec. IIVBJ) conjugation symmetry for o~(k), and subsequent results for this 
quantity apply only for k <C N . By the usual rules of Fourier transforms this replacement of a sum by an integral 
requires us to turn the sum over k of the inverse transform in Eq. @ into a sum with infinitely many terms. This, 
as we will see from explicit GO results, leads to dropping non-scaling terms in p(x) which are exponentially-small in 
N and is harmless. 

Then we go further and replace the infinite sum over k by an integral over k. We will see in sec . IIV 51 that for any 
given form of cr(k) we must modify p{x) to a new function p(x), while at the same time changing the integral over 
x that expresses o~(k) by an integral with infinite limits. The modified function p(x) does not obey p{x) = p(l — x), 
but is still even in x, which is enough of conjugation symmetry for our purposes. We see by studying GO's explicit 
results that p{x) is a certain limiting form of p(x) 7 valid when the parameter 



is small compared to unity. We show in sec. HVCl that the smallness of e is a form of the DGA that is useful, since it is 
compatible with the scaling oip(L) as given in Eq. Ijljl. It turns out that requiring e <C 1 is appropriate for justifying 
the use of Fourier integrals rather than sums. Taking the usual fundamental string tension, and A -1 ~ M ~ 600 
MeV, we find that e is about 0.08. 

Having made the kinematic simplifications alluded to above, we explore crudely the behavior of the scaling density 
p(x) in the limits of small and large x. This behavior is compatible with, but does not necessarily imply, a linearly- 
growing string tension. In fact, the behavior at large flux is precisely that found in |llj | using the DGA, but our result 
is compatible with scaling [Eq. 10}] . 

Sec. IVII summarizes the paper, and in addition there are two appendices covering details of center- vortex actions. 

In future works we will make more extensive studies, both numerical and analytic, of the master rate equation 
that determines the p(L). But even at this primitive stage of development, we hope to convince the reader that it is 
plausible that the center-vortex picture leads to <r(k) ~ fc for small k. 

II. DESCRIPTION OF CENTER VORTICES AND NEXUSES IN THE GAUGE THEORY 

Throughout this paper we use gauge potentials that are anti-Hermitean matrices constructed from the canonical 
potentials and multiplied by the gauge coupling g. As discussed in the introduction, we write explicit formulas in 
d — 3 Euclidean space, but the transcription to d = 4 is straightforward ; it leads to an unconventional string theory 
rather than a field theory. 

Center vortices are (quantum) solitons of the effective action, and they are essentially Abelian objects that can be 
superposed. An elementary center vortex is described by specifying a closed curve and a group matrix. The group 
matrices must be chosen so that the magnetic flux F, defined through the Wilson loop as 

— -Tr fc exp[<f dziAi] = exp[ikF], (7) 

is an element of the center. Here Tr k indicates a trace in the k representation, of dimension D k - (Of course, k = 1 is 
the fundamental representation.) If the Wilson loop T is linked exactly once in a positive sense with the vortex curve, 
then F is the flux associated with the vortex. It must have the form F = 2ttJ/N, where J is an integer from 1 to N, 
in which case we speak of a J vortex. 

The idea behind our picture [2j of center vortices is that infrared slavery requires the generation of a dynamical 
gluon mass |2j]. This, as seen in Schwinger-Dyson equations for the gluons, induces a gauge-invariant mass term that 
is simply a gauged nonlinear sigma model. The effective action for the gauge theory is then 



1 /■ „ , 1. 



l eff 



d 3 x[--TrGf 1 -M 2 Tr(U- 1 D l U) 2 } (8) 



2 / ~ ~L q ' "13 



g 2 J L 2 



where U is an N x N matrix in the fundamental representation, transforming as U — ► VU when Aj — ► VAiV^ 1 + 
VdiV~ l , and g 2 is the square of the coupling constant. This matrix U is to be eliminated by minimizing the action in 
which it appears. If the mass M were considered as a fundamental object, surviving unchanged into the ultraviolet, 
this effective action would not be renormalizable. In fact, one knows |22( that the mass is a function of momentum q 
decreasing at large q 2 as 

const (G 2 A 

M 2 (q) - ^. (9) 

q- 

Some such decrease is required for the Schwinger-Dyson equations to be solvable. We will only use the effective action 
in the infrared, so this complication of the ultraviolet behavior of the mass term is irrelevant for us. There is one 
exception: if the mass M is wrongly taken to be constant, the mass term in I e ff has a logarithmic divergence. Of 
course, this is not real, and we will ignore it. 

An elementary unit-flux center vortex, a solution to the equations of motion of the effective action, has a form 
essentially equivalent to a Nielsen-Olesen vortex: 

4"] (a,) = ^Le ijk dj I dz t [A M (z - x) - A a (x - z)]. (10) 

Here Q r is one of a set of N matrices defined below, whose properties are detailed in Appendix A, and Am(o) is 
the free Euclidean propagator of mass M (0). Provided that the curve 7 is closed, the Aq term is pure gauge, and 




FIG. 1: Center vortices merging and dissociating. A J = 1 vortex (labeled 1 for Qi) meets a J = 2 vortex (Q2 + Q3) at the 
bottom vertex to form a J = 3 vortex, which later dissociates at the top vertex. 



corresponds to the contribution of the U matrix. (If the curve were not closed, long-range monopoles would sit at the 
ends of the curve.) The field strength of the vortex decreases exponentially away from the vortex curve 7. 
The matrices Q r are 



diag — , — , . . . — , 

5V iV' N' N' 



1 1 

~ 1+ iV' N' 



1,2, 



,N 



with the -1 in the rth position. Each of these matrices obeys 



exp[2iriQ r ] = exp[27ri/iV] 



(11) 



(12) 



and so has unit flux. Since each Q r can be transformed into any other by an element of the permutation group, one 
may think of the index r as a label for group collective coordinates of the vortex. These matrices have the property 
(see Appendix A) that the sum of all N of them vanishes. 

One may immediately generalize this by adding any number of solitons of the form of Eq. I|10[) . In particular, if for 
a given curve 7 we use the flux matrix 



F {J} - Qh 



■ Qi 



(13) 



with all Q indices distinct, we find that the flux of Frjy is J. We can construct a conjugate vortex with flux matrix 
F{n-j} by adding together the N — J distinct Q matrices with indices different from all those in Frn of Eq. i|13|l . 
Because the sum of the Q's is zero, it follows that 



F{J} + F {N-J} — 0, F[_jy = -F^jy - F{ N _jy. 



(14) 



So antivortices of flux —J are entirely equivalent to vortices of flux N—J, illustrating the mod N nature of flux addition. 
As before, the matrix indices on Fijx can be shuffled by permutations, which yield group collective coordinates (see 
Appendix A) . A vortex of flux J is identical in form to that of the unit vortex in Eq. H10|) except that Qj is replaced 
by F{jy. In what follows we suppress the collective coordinate index [r] in Eq. (jlOjl and similar equations. 

Note that whatever the flux J is, vortex strings always have the same thickness of M _1 . 

It is important for the physical model discussed later to recognize that center vortices can merge and dissociate. 
We describe this by a vortex soliton of the form 



-4, 



2tt 



tijkdj 22 



dz(a) k F{j a }{A M [z(a) - x] - A [z(a) - x}} 



(15) 



where the closed curves labeled by a can have segments in common, as illustrated in Fig. ^ 

Note that flux is conserved in processes of merging and dissociation as shown in Fig. ^ There is another type of 
merging allowed, in which flux is only conserved mod N. We begin by writing a partial expression for the associated 
soliton: 



2tt N f 

Ai(x) = -r€ijkdj'Y] / dz(a) k Qa{A 

1 a=l J 



M [z{a) -x)- A [z(a) - x]} 



(16) 




FIG. 2: Center vortices of total flux J — N can be annihilated or created from zero flux. The labels 1, 2, ... correspond to the 
labels of Q matrices. 




FIG. 3: A recombination event, in which vortex 2 (that is, of flux matrix Q2) combines with vortex {13} at the bottom, and 
dissociates into vortex 1 plus vortex {23} at the top. Black circles are nexuses or antinexuses which change 1 into 2 or vice 
versa. 

where the sum now is over curves a which all begin at the origin. Because the sum of all Q a is zero, there are no 
long-range monopole fields associated with this termination of the curves. 

It is still necessary to close the curves somewhere away from the origin. One way is simply to repeat the process in 
Eq. (|16fl at some other point x — u. This is illustrated graphically in Fig. [2] Note that this is topologically precisely 
the same as a bayonic Wilson loop in SU{N). 

There is one other major type of configuration remaining; it involves nexuses 20] . An elementary nexus is a point- 
like region on an elementary (unit flux) vortex string where the flux matrix changes smoothly (on the scale set by the 
mass M) from Qi to Qj. It is truly a non-Abelian object, with an explicit description too complex for us to record 
here; it will be enough to know that the nexus is essentially a point on a vortex string which changes the Qi around. 
A nexus (or antinexus) conserves the flux J mod N but changes the matrices Qi which make up this flux. The flux 
difference Qi — Qj is essentially a root vector of SU(N); it has one diagonal element of 1, one element of -1, and 
the rest are 0. This is the Pauli matrix 03 for some embedded SU{2) of SU{N), and constitutes the flux matrix for 
what is basically a 't Hooft-Polyakov monopole whose field lines are deformed into the strings of the center vortices 
on either side of the nexus. Since every vortex string is closed, every nexus must be accompanied by an antinexus 
somewhere else on the vortex string. In d = 4 the linking of nexus world lines to vortex surfaces leads to topological 
charge occurring in localized lumps of charge 1/iV, but globally of integral total charge. This is not of concern to us in 
three dimensions. But nexuses do lead to a simple description of important recombination processes, a little different 
from the simple merging and dissociation of Fig. ^ An example is shown in Fig. |3 showing a Qi vortex exchanging 
with a Qi vortex during merging and dissociation. 



III. 



A THREE-DIMENSIONAL SCALAR FIELD MODEL OF NEXUSES AND CENTER VORTICES 



It has long been known |8| that the partition function of simple closed strings in d — 3 can be described by 
the partition function of a scalar field theory. The essential ingredient is that the Feynman-Schwinger proper-time 
representation of the (trace of the logarithm of the) free scalar propagator has an interpretation in the language of 



closed strings: 

, rA . f°° ds f , , . . f s , .Mi . 9 M 2n , , , 

Trln[A M ]=J — f(dz)exp{- J dr\-fz % + -f]}. (17) 

As a stringy object, in the above integral the parameters s,t are not proper times but rather physical lengths along 
the string, and (dz) is an integral over all closed string paths. The masses M 12 are mass parameters derived from 
some underlying theory; they are related to the physical mass by M 2 = MiM 2 . 

For oriented strings, the string partition function is an integral over all closed random walks and lengths: 

^r ing =j:^{f^f(dz)eM-[dr(^ + ^)}} N . (18) 

Here the z 2 term is the familiar action for random walks of length s. It involves a mass factor Mi related to the 
diffusion rate of the random walks. The other mass M 2 gets positive contributions from string self-energy per unit 
length and negative contributions from entropy, so it is possible that M 2 and M 2 are negative, and the field is 
tachyonic, signaling the formation of a string condensate. In that case, the free theory does not make sense, but 
stabilizing interactions, such as a four-point term in the field theory, must be added. Then the scalar fields pick up 
vacuum expectation values, and their masses are physical. It turns out that our scalar field theory has no continuous 
global symmetry, so there are no Goldstone bosons. 

The factor 1/s in the "proper-time" integral corrects for overcounting strings with different starting points on the 
closed curves. For unoriented strings, this factor is l/(2s). The sum over N gives the usual free-field scalar partition 
function for oriented strings: 

% eW = exp{Trln[A M ]} (19) 

To represent the particulars of the gauge solitons of sec. [n]we introduce N scalar fields <pi, i — 1, 2, . . . N. There 
are relations among these fields, following from the properties of the Q matrices. The first relation is conjugation: 



?N-i 



(20) 



For N even, the field </>jv/ 2 is self-conjugate. The other relation is that one of the fields, say <pN, can be expressed in 
terms of the others, which corresponds to the vanishing of the sum of the N Q matrices. In the scalar field language, 
this amounts to writing a vertex in the action of the form 



iV 



G d 3 xY[<t> t (21) 



that corresponds to the vertices in Fig. @. 

It is possible in principle to calculate the properties of the scalar-field action directly from the gauge theory. For 
example, the contribution of one unit-flux vortex soliton, as displayed in Eq. I |10| l. to the free scalar-field action Jo is 
found by calculating the gauge action of Eq. © based on the solitonic potential of Eq. (|10|) . We record here only 
the B 2 part; the mass term gives a similar term. 

M 3 7T f f 

h = —Tr{Q 2 ) dr dr'zi(T)zi(T')ex.p{-M\z(T) - z(t')\}. (22) 

(We do not need to indicate an index on Q; all Qs give the same result.) The exponential term comes from the 
integral over all space of the square of the magnetic field strength Bi, which is [up to constants already absorbed in 
Eq. (|221)] 

Bi(x) = f dz l M 2 A M {z - x). (23) 

In the infrared regime, where all length scales are large compared to M , one can make contact with canonical 
field theory by making the large-M approximations 

z t {r)z t {r') ~ z t (r) 2 , \z{r) - z(r')\ ~ \z(r)\\r - r'\. (24) 

The integral over r' can now be done, and (omitting exponentially small terms) results in 

h=Tr{Q 2 )^- J dTe~ l (T)z t {r) 2 (25) 
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where the einbein e is just |£(t)|. This contribution to the free action is reparametrization invariant; the form quoted 
in Eq. H17|l is the special case |i(r)| = 1. One can add to this kinetic action a term representing energy or entropy 
per unit length, of the reparametrization-invariant form J dTe{f)M%]2. By comparison of Eq. Q17|l and Eq. 125|) one 
sees that the mass Mi is 2irTr(Q 2 )M 2 /g 2 . As is well-known, in d = 3 all dynamical masses are linear in g 2 , and 
therefore so is M\. 

One must be careful to understand that the large- M or infrared form of the action in Eq. I|24(l is not appropriate 
to study the ultraviolet behavior of the scalar theory. The contribution of one closed vortex, of whatever shape or 
length, to the free scalar propagator A Qc f, is the integral of exp[— 1 — (M 2 s/2)] over all curves and lengths: 

p°° r x ( m r s f s ~i 

Ao*(a;)=/ dse- Ahs / 2 (<fe)exp-{-^/ dr d T / i i (T)ii(T , )ex P [-M|«(T) - z(r')\]j (26) 

before any approximations are made. By writing the integrals over r and r'in the exponent as the limit of a sum of 
terms one can check that no terms singular at s = arise from the path integrals. Normally, the path integrals lead 
to a factor s~ 3 ' 2 in the s integral, which is responsible for ultraviolet divergences in the propagator. But no such 
terms arise in our case, because the scalar field theory is built on solitons, rather than on pointlike fields. The scalar 
propagator behaves at least as well as M 2 k~ 4 for large k, which means that even the iV-point vertices discussed below 
lead to convergent Feynman graphs. 

The non-local (in proper time) corrections following from Eq. (|26H to the standard local propagator of Eq. 117|) 
can be expressed as a power series in the deviation of the trajectory z(t) from a straight line. These are not easy 
to characterize, because even for a simple circular orbit the non-local terms lead to a divergent series that is not 
Borel-summable. If we specify the orbit as part of a circle 

z(t) = %os( J), sin(^), 0] (27) 

the double proper-time integral becomes 

J dr /dT'cos(^-^)exp{ -2MRsin[ * T ~ * ]}. (28) 

The expansion of this integral in powers of 1/(MR) is a factorially-divergent series with all terms of one sign. This 
expansion is closely related to the asymptotic expansion of the Bessel function Iq(2MR). We will not explore this 
subject further in this paper. 

We next discuss the "baryonic" vertex of Fig. @ and Eq. (|21|) . Let us calculate the contribution Ig to the 
logarithm of the scalar partition function Z^ coming from the solitonic configuration of Fig. |2J) [see Eq. (|16j . in the 
infrared or local approximation. We assume the strings in Fig. J3J are well-separated on the scale of M , except for 
the two regions of size M _1 where they all meet. Using the same techniques as for the propagator of a single scalar 
field we find 

In Z - e- /N J d 3 x J d 3 y[A 4x - y)] N + ... (29) 

where In is the contribution to the action coming from the regions of size ~ Af _1 near the points x, y where all N 
strings come together. The factor exp[— In] is a factor of G 2 , where G is the iV-coupling constant in Eq. J2U- 

Evidently the Feynman graph represented by Eq. (|29|l is exactly the same as the graph of the solitonic strings 
in Fig. J5J. And this is true generally: Every graph of solitonic strings has a counterpart Feynman graph in the 
scalar field theory. The sum of all connected Feynman graphs represents the logarithm of the solitonic- string partition 
function. 

In particular, nexuses amount to a quadratic term in the action. The most general quadratic term is 

i 3 xJ24v^ r (30) 

The diagonal terms represent single-string energy and entropy per unit length, which is the same for all vortices, and 
the off-diagonal terms represent the nexus-mediated flipping of Qt to Qj . The nexus action is the same whatever the 
values oii,j. As a result, the matrix V has one element rr? d on the diagonal and one element m^ for every off-diagonal 
entry. Such a matrix is easily diagonalized, and has N — 1 eigenvalues m 2 — m 2 w and one eigenvalue of m 2 , + (N— l)m 2 . 
The N — 1 eigenvectors of the first eigenvalue are of the form J2 x i4>i with J2 x i = 0, and for the second eigenvalue 
the eigenvector is ~Y^4>i- In the large- N limit this eigenvector does not propagate, because its eigenvalue is of order 
N. This is consistent with the fact that in the gauge theory there are only N — 1 unit-flux solitons. 



Finally, it is worth noting that the complete lack of ultraviolet divergences, plus the fact that all masses are linear 
in the gauge coupling g 2 , is all that is needed to prove the exact result of [23|], which is that the full quantum effective 
action T, as a function of the zero-momentum condensate (G 2 ■), has a minimum corresponding to a non-zero value 
for this condensate. 

Before going on to the rate equation derivable from this scalar field theory, we first discuss the relation of the 
Wilson-loop VEV to probability densities for vortices piercing the loop. 

IV. WILSON LOOP EXPECTATION VALUE AND AREAL DENSITIES 

We use the the GO formulation of Greensite and Olejnik [12], with some minor additions. This formulation goes 
beyond the usual DGA formula for relating a(k) to the discrete Fourier transform of the areal densities and avoids 
the problem, mentioned both by [ll| and GO, that if the standard DGA formula is used to relate a string tension 
linear in k to the vortex areal densities, these densities are singular and do not show correct large-N behavior. We 
will find a modified approach to the DGA based on GO's work which avoids this difficulty. 

A. Areal densities and the Wilson loop VEV 

Consider a large planar Wilson loop of area A, divided into many squares whose size is given by the vortex correlation 
length A. We expect A ~ 1/M. Let p(L) be the probability that a vortex of flux L exists on any square, so that the 
areal density of vortices is p(L)/X 2 . We assume that only one vortex can sit on any one square, since (as discussed 
in the next section) two vortices attempting to occupy the same square will either repel each other or attract each 
other to form a new single composite vortex. On the average there are as many antivortices as vortices, and so the 
antivortex probability p(—L) = p(N — L) is precisely equal to p(L). (So on the average the net flux through any 
Wilson loop is zero. The area law of confinement arises from fluctuations in the total vortex flux.) 

The Wilson loop contains B = A/A 2 squares, and we calculate the VEV in terms of the probabilities in the k 
representation as 

/2itk 
dxiAi]) = P (0) B +Bp(0) B - 1 2p(l)cos( — )+ (31) 

+ B^-l) p(o)B _ 2[2Mi)cos( 2^ )]2 + SK0) b- %(2)cos( 2^ ) + 

where D^. is the dimension of the representation k, and we have used the equality of vortex and antivortex probabilities. 
The presence of the term B{B — l)/2 rather than B 2 multiplying p(l) 2 says that if two vortices are present they must 
occupy different sites, for reasons given above. 
The sum in Eq. (|3T|) is 

{Y / p(L)exp[^^}} B = {p(0) + 2Y / P(L)co S (^ r )} B (32) 

L=0 L=l 

where the prime on the sum indicates that it goes from 1 to (N — l)/2 if N is odd. If N is even, the sum goes 
to (N/2) — 1 and one must also add p(N/2) cos(7rfc) to the sum. This sum is supposed to equal the area-law result 
exp[— a(k)A], where a(k) is the string tension in the k representation. Since B = A/A 2 this leads to, as in GO and 
Eq. © 

/ N — l 

expha(fc)A 2 ] = p (0) + 2^p(L)co S (^)= ^p(L)exp[^p]. (33) 



L=l 



The behavior of <j(k) is constrained by conjugation symmetry [p(L) = p(N — L) = p(—L)]. For k restricted to integer 
values we find 

a(k)=a(-k)=a(N-k). (34) 

This furnishes an extension of er(fc) to negative integers. A simple change of variables, using conjugation symmetry, 
shows that Eq. l|55|l can be written as 

N/2 2-nkT 

exp[- ( x(fc)A 2 ] = Y, P(i)^[— ]• (35) 

-N/2 
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[Here the limits ±JV/2 are nominal but apply for large N; the real limits can be worked out from the definition of the 
primed sum in Eq. (|32|) .] 

B. The large-N limit and conjugation symmetry 

We will be interested in the large- N limit, which has certain subtleties that we explore next. 
The inverse transform corresponding to Eq. (|33l is [Eq. Q]: 

P^ = J? E ex P[^^l exp[-a(fc)A 2 ]. (36) 

fc=0 

As long as L is an integer it is straightforward, using conjugation symmetry as in Eq. (|34|l . to show that the inverse 
transform can also be written 

PW = ^ E exph^WA 2 ] cos[^]. (37) 

-N/2 

[As before, the limits on the sum are nominal.] From either of these equations it is very reasonable to expect that the 
p(L) have the large-A scaling property [already noted in the Introduction] 

with conjugation symmetry implying p(x) = p(l — x) — p(—x). If so, then it is a good approximation to replace the 
sum in Eq. I|33|l by an integral: 



exp[-(j(fc)A 2 ] = f , 
Jo 



dxp{x)e^ lkx (39) 

The error made is at least as small as 0(1/N), and in certain cases of physical interest is exponentially-small in N. 
Note that we have lost conjugation symmetry for a(k), but the property er(fc) = cr(—k) which follows from Eq. I|39|l 
and conjugation symmetry for p{x) replaces it for all practical purposes. As long as k is a positive or negative integer, 
conjugation symmetry for p(x) allows us to write this equation in the alternative form: 

,1/2 

exp[-cr(fc)A 2 ] = / dxp(x) cos[2irkx\. (40) 

7-1/2 

At this point we will use Eq. 140|) to define cr(k) for general non-integral k, when we need to do so. 

So far we have replaced one of the discrete Fourier transforms, Eq. I|35(l . by an integral over a finite domain, 
Eq. I|40l) . This means, of course, that the inverse Fourier transform, Eq. i|3t)|) . remains a sum over integral values of 
k, but the sum extends to k — ±oo: 

oo 

p(x) — 2, ex P[ — "Zirikx] exp[— cr(fc)A 2 ]. (41) 

k— — oo 

One might now ask whether it is legitimate to replace the sum in Eq. I|37|l by an integral with limits ±oo, provided 
that the sum over k is well-behaved. This integral form would explicitly exhibit scaling symmetry. But Eq. I|40|l and 
this integral would not, in general, form a Fourier transform pair. This means that for a given a{k) the result of 
replacing the sum in Eq. (|37|l by an integral defines not, as one might hope, p(x), but another function p(x), which 
we hope is closely related: 

/>oo 

p(x) = / dke- a{k)x \- 2 * lkx . (42) 



The general principles of Fourier integrals tell us that the true Fourier inverse that expresses exp[— <r(fc)A 2 ] is found 
from replacing the limits ±1/2 in Eq. (|40|l by ±oo and p(x) by p(x): 

/oo 
dxp(x) cos[2nkx]. (43) 

-oo 
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Now both cr(fc) andp(cc) are defined for certain complex domains of their arguments, and both are even functions. The 
original conjugation symmetry for both of these functions has been lost; there is no reason, for example, to believe 
that p(x) = p(l — x). But it is manifest that p(x) — p(—x). 

Is there much of a relation between p(x) and p{x)l The answer is that there can be, if in Eq. 1|4(J[1 the integration 
over x receives dominant contributions only from \x\ <Cl/2, in which case ±1/2 is almost as good as ±00 [the new 
limits in Eq. (|43(l . As we will now see, the explicit formulas of GO furnish an illustration of this point, in which 
smallness of the parameter e introduced in Eq. © is equivalent to effective domination of the integrals in Eqs. (|4(JI43|) 
by small values of x. We will also see that smallness of e is a form of the DGA. 

The explicit phenomenological example of GO is 

N N 

cr(fc) = fcer(l), k < — ; =(N- fc)cr(l), k> —. (44) 

They then calculate the corresponding p(L), which has terms exponentially-small at large N. Such terms are auto- 
matically dropped by using the discrete Fourier transform of Eq. (|41|l with infinite limits, and we find 

1 - 7 2 
P ^ = (1 - 7 ) 2 + 2 7 [1 - cos(27rfca;)] (45) 

where 

1 = e-< y ^ x2 =e- 27re . (46) 

Consider now what happens to GO's p(x) [Eq. 1|45|) ] when both e and \x\ are small. The result we will identify with 
p(x): 

7r(e z + X ) 
The probability p{L) corresponding to p above is 

p(i) = AN^TlXy (48) 

This p(x), used in the Fourier integral of Eq. I|43|l . yields a string tension 

exp[-er(/c)A 2 ] = cxp[-27re|A:|] (49) 

which is the GO input string tension [Eq. i|44[l ]. provided that k < N/2. And if the string tension is defined by 
Eq. H49[) for all k, the Fourier integral of Eq. 1)43|) yields, as it must, p{x). The importance of small e is that the poles 
of the integrand of Eq. I)43|) are at x = ±ie and thus far from ±1/2. 

C. The DGA 

Let us close this section by relating the smallness of e to the DGA. GO have observed that scaling makes it very 
difficult to understand how the DGA could be applicable, since p(0), the probability of no vortex in a A square, is of 
order 1/A^ and so there must be vortices almost everywhere. In the formal large- N sense this is certainly true, but for 
fixed N, however large, it is always possible to choose e > l/(xN) so that p(0) from Eq. (|47|l is close to unity while 
the next-largest probability, which is p(l), is an order of magnitude smaller. [One might also note that the formal 
e — > limit of p(x) in Eq. I|47|) is S(x).] Of course, in the real world of physics we cannot choose e, and in any case 
we would not expect the real world to be describable in DGA terms. Nevertheless, it is possible to envisage a DGA 
that both accommodates scaling and allows for the areal densities of vortices to be fairly small. 

The usual form of the DGA comes from subtracting unity from each side of the basic GO result in Eq. (|33fl . and 
using the normalization of probabilities to eliminate p(0): 

1 - exphcx^A 2 ] = 2 5>(L)[1 - cos(-^-)]. (50) 

L=l 

At least formally, the DGA follows from assuming er(£;)A 2 is small and from saving only the linear term in this quantity 
on the left-hand side of Eq. I|5U|) . Since in our example a{k)\ 2 = 27re|fc|, this would require 27re (and 2nke) to be 
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small. As mentioned in the introduction, ref. [Tj uses this conventional DGA to study what probabilities p(L) are 
needed to reproduce various forms of er(fc). If tr(fc) ~ fc for N s> k, Ref. 11] finds [in effect by linearizing the left-hand 
side of Eq. (SJ)] that for L > 1 



N 

I 2 



p(L) ~ 75 . (51) 



But, as [llj and [T^l note, it is hard to understand this formula for general L > 1, because an areal density proportional 
to A~ at small (i. e., non-scaling) values of L would indicate a pileup of many vortices on a single square. The simple 
change from the denominator L 2 of Eq. (|51|) to L 2 + N 2 e 2 in our Eq. (|48|l yields a well-behaved scaling form for p(L) 
at all L, including 0, while maintaining the behavior L~ 2 in the regime eN <I< N/2. 

The basic point for us is that the simple form of p(L) in Eq. IJ4"7|) (1) scales as it should; (2) leads to p(Q) ~ 1/N; 
and (3) has an asymptotic behavior L~ 2 which is compatible with a fc-string tension rising linearly in k for small 
enough k. This form was discovered by replacing the limit N or N/2 on certain flux sums by oo, and it is just this 
replacement we will make below in our approximate rate equation. It will then be seen that this rate equation has a 
solution with all three properties listed above. This is not to say that the solution of the rate equation is exactly as 
shown in Eq. (|47|l . as discussed later on. But our present limited investigation of the rate equation reveals no obvious 
imcompatibilities with Eq. 147|) . 

V. THE MASTER EQUATION 

The master equation for vortex dynamics describes the evolution of vortex areal densities as a function of a variable 
z, which can be thought of as giving the intercept of a test plane (essentially a planar Wilson loop), perpendicular to 
the z axis, with that axis. For a given configuration of branched vortices in three-space, as z varies different cross- 
sections of the configuration are seen (see Fig. 0J. These cross-sections look like a gas of d = 2 molecules merging, 
dissociating, and recombining; the molecules are made of atoms that are the cross-sections of unit-flux vortices, and 
they attract and repel according to laws given below. All such processes strictly conserve flux number L mod N . A 
flux matrix -F{l} is a sum of distinct Qi matrices [(see Eq. (|13|) ]. so we may think of any vortex of flux L as made of 
constituents, which we call atoms, that have unit flux and are labeled by the indices on the Qi. A center vortex of 
flux L ^ 1 is a molecule, made of those atoms which appear in the flux matrix FrL}. We can always take flux values 
to lie in the range 1 < L < N — 1, and vortices with L > N/2 are antivortices. On the average, as many vortices of 
flux L as of flux N — L pierce any sufficiently large plane, so the net flux is equivalent to zero. We assume ergodic 
behavior, so that averaging over z is the same as averaging over many different vortex-nexus configurations. 

The atoms' and molecules' motion is essentially diffusive, because they are cross-sections of random walks. Every 
center vortex, of whatever flux, has the same size, roughly equal to the correlation length A of the vortices, so we can 
imagine the test plane divided into squares of area A 2 . Each square can hold one vortex of any flux L, with probability 
p(L) as described above. The probability p(0) is the probability that there is no vortex at all in a square. 

We are interested in large N, first because in any case N must be greater than three or there is no fc-string problem, 
and second because we would like to be able to enter the scaling regime, defined in Eq. (|38|l . that contains vortices 
whose flux L scales with N as N grows. As we have already mentioned, even with scaling it is possible to adjust 
parameters so that the probability of no vortex, p(0), is considerably greater than the probability p(l) of one vortex. 
The scaling regime is separated from this regime of sharply-varying p(L) at a flux L ~ eN , where e <?C 1 is (t(1)A 2 /(27t), 
as before. The discussion in the previous section tells us that replacing flux sums by integrals and taking certain 
limits of these integrals to be infinite is qualitatively appropriate, provided that e is fairly small. 

We will use the flux dependence of the semiclassical actions of various vortices to decide whether processes involving 
an action change is allowed or not. That we can use the semiclassical actions of center vortices at large N is certainly 
an assumption which can be challenged, on the grounds that the actions in question diverge as N gets large [see 
Eqs. H22|). HA12fl ]. But one knows [24J that the group entropy for center vortices can completely cancel the action 
at leading order in N, and that if this happens, the leading large- TV corrections preserve the functional form of the 
action except that this action is multiplied by const /N 2 . The resulting then scales at large N. However, the fate of 
non-leading large- N corrections has not been studied, and we have nothing to say about them. In any case, we must 
and do assume that some such phenomenon of cancellation of action and entropy takes place, or else center vortices 
do not survive the large- N limit at all |25|. 

In any case, we make the assumption that the probabilities for particular merger and dissociation processes are 
large or small according to the semiclassical action difference associated with the initial and final states of merger 
or dissociation. One might expect, as in the chemical law of mass action for rates, that rates depend on exp[— AJ] 
where A/ is a difference of actions. The rationale for this can be seen in, for example, Eq. (|29|l . where the "baryon" 
vertex coupling G is of the form exp[— /]. But the law of mass action strictly holds only for dilute systems, and we 
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FIG. 4: A test plane perpendicular to the z axis is moved along that axis, revealing a continuous picture of merging and 
dissociating vortices piercing that plane. Shown are two vortices recombining and another vortex dissociating as seen with z 
as the "time" . 

take a different approach. In fact, we will drastically simplify the rate equation by allowing flux-conserving processes 
for which the action difference is favorable, and forbidding those for which it is not. The actions involved depend on 
quantities such as Tr(F{jyF{xy), and comparing these actions for initial and final states of vortex processes yields a 
set of rules for probabilities of these processes, plus force laws which predict whether particular vortices will attract 
or repel each other. Next we summarize allowed and forbidden processes in a set of rules for action and for forces 
between vortices. 



Force and action rules 



A matter of notation: From now on, when we refer to an action, we mean (unless the context makes clear otherwise) 
a trace quadratic in flux matrices, without factors of AirM/g 2 , etc. So wc define 



I 



{£} 



TrF ?L}- 



(52) 



Given the assumption that it is meaningful to compare the semiclassical actions, and the traces (see Appendix B) in 
question, we will now see that center vortices of flux L can be understood as a composite of L unit-flux vortices which 
have certain laws of attraction and repulsion. Generally, a vortex of flux matrix Qi is repelled by another vortex of 
the same type, and attracted to vortices of different type but by an attraction smaller by a factor of 1/N. (This is to 
be expected for some sort of rough force balance, since there are roughly N times as many ways for a vortex to see 
another of different type as ways to see one of its own type.) 

In merger processes it is possible to form flux matrices where some Qi appear more than once. If some indices are 
repeated, we say that there is an overlap. As in Appendix B and Eq. Q13p. we use the notation F{ L y for the sum of 
L matrices Qi with no indices repeated. If R of the indices are each repeated once and only once, we indicate this 
by -F{z, ; _r}. Of course, this is an incomplete notation, since in general the precise content of the set {R} of repeated 
indices must be specified. Furthermore, one could imagine that some indices are repeated once, others repeated twice, 
etc. We restrict ourselves to single repeats because the more Q matrices of like type in a flux matrix, the higher the 
action and the less likely this is to occur. 

The actions I{l,r} calculated as a trace in Appendix A, Eq. (|A13|I . are exactly reproduced by the following rules 
for calculating the self and mutual actions of unit-flux constituents: 

1. Each constituent (unit-flux vortex) has self-action (N — 1)/N. 

2. Unlike constituents have action —2/N per interacting pair. 

3. Like constituents have action +2 per pair. 
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Such rules of attraction and repulsion automatically lead to a stable system, since every action in the system is the 
trace of a real matrix squared. 

As an example, the action J/j,\ of a vortex of flux L is L(N — L)/N if there is no overlap. We interpret it as 
follows. Each of the L constituents has self-action (TV — l)/N, and is attracted to unlike constituents with an action 
contribution —2/N. As a result, the action of a vortex with no overlap is 



L( ^ + V-i)rf) = ^= x) 

N 2 V A N ' N 



hL} = ^^ + *L(L - 1)(^) = ^V^- (53) 



If there is overlap R, one must add 2R to this action, according to Appendix B, Eq. l)A13jl . This positive addition 
indicates that like constituents repel. 

We can transcribe the action and action difference associated with various processes into rules for attractive and 
repulsive forces between vortex atoms. Consider the sum of two vortices in the form given in Eq. l|15f) . and suppose 
that the string for vortex a = 1 is a straight line along the z axis. The string for vortex a = 2 is the same, except 
that it is displaced in the xy plane by a distance r. Up to factors irrelevant for the present discussion, the action of 
the pair is easily found [using Eq. (jAllfl ] to have the form 

/ = I {Kl} + I {K2} + 2e- Mr [-^l + R] (54) 

where R is the overlap between {i^i} and {-ft^}- Evidently, if the overlap is small enough the third term on the right 
is attractive, but otherwise it is repulsive. We will use this force rule and the action rules in the master equation. 

B. Qualitative aspects of the master equation and its solution 

As discussed above, the master equation is a rate equation in the variable z, where z is the point on the z axis 
where a test plane of fixed contour and area A, perpendicular to this axis, intersects it (see Fig. |3J|- 

In general, there are a large number of many-vortex interactions. To keep things manageable, we concentrate on a 
limited number of these, for which we furnish simple go/no go rules. They are: 

1. Merger of two vortices into one 

2. Recombination of two vortices into two other vortices 

3. Decay of one vortex into two vortices 

The rules for these processes are based on the action and force rules of the previous subsection. In implementing the 
rules we will replace sums by integrals and set some upper limits in these integrals to oo rather than to N, in much 
the spirit of sec. II V Bl 

The most important action or force rule is that two vortices with zero overlap tend to attract, and that the state of 
lowest action is the single vortex which combines all the constituents. Furthermore, if they have finite overlap but still 
attract, they will, in whatever recombination of constituents occurs, prefer to separate the overlapped constituents 
entirely. That is, suppose vortices with {J} = {123456} and {K} = {12789} merge to form the virtual vortex 
{J + K;R} with {R} = {12}. Lower-action states are available, so this virtual vortex will decay quickly. The decay 
state that conserves constituents and has the lowest action is the sum of vortices {123456789} and {12}. 

Consider now a merger process in the scaling regime, in which vortices with fluxes {.Ki} and {K2} meet at a point 
(that is, they occupy the same correlation- length square in the test plane). Neither vortex has any internal overlap 
(otherwise it would dissociate, as said above), but when merged there is an overlap. The probability distribution 
of overlap is described in Appendix B. The average overlap of the configuration is (R) = K1K2/N, and it is in the 
scaling regime if Ki, K^ are. The merged configuration {K\ + K^\ R} has action [in the sense of Eq. (|52|l ] 

KiK 2 , , . 

I{k 1+ k 2 -b,} = I{i< 1 } + I{K 2 } + 2[R —J. (55) 

On the average, the overlap repulsion just cancels the attraction between the two vortices, so the merger is action and 
force neutral. However, by recombining constituents, the action can be lowered. The general rule is that the most 
favorable state is a single vortex with the combined flux of the two original vortices, as one sees from Eq. <|A12|1 . 
However, if there is overlap, one must arrange it so that Qi of the same type go to different vortices. 
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If the initial configuration is rearranged into {L} + {K\ + K 2 — L} in such a way that neither {L} nor its partner has 
internal overlap, the action of this recombined configuration, once the two vortices have separated several correlation 
lengths, is 

I{L} + I{k 1+ k 2 -l } = I{K l} + I{K 2 } -Jj( L - K i)( L - &*)• (56) 

The go/no go rule is that if the product [L— Kx)(L— K 2 ) multiplying 2/N in this equation is positive, the recombination 
process is allowed, otherwise disallowed. One easily sees that the conditions favoring the recombination are 

L < min(Ki, K 2 ) or L > max(K 1 ,K 2 ). (57) 

In addition, for a process to be allowed we require that there be no overlap in either the {L} vortex or the {K\+K 2 — L} 
vortex, so that 

L>X^,K 1 + K a -L>*g>. (58) 

It is important to note that the constraints of Eqs. I|57l I58JI are conjugation-symmetric. Conjugation symmetry 
requires invariance under the exchanges 

L^N-L, Ki-tN-Ki, K 2 ^N-K 2 . (59) 

Clearly, the conditions of Eq. (|57|) are invariant under conjugation symmetry, since they derive from the conjugation- 
invariant combination (L — K\){L — K 2 ) of Eq. (|56() . And we note that the second inequality of Eq. I|58|) . stemming 
from no overlap, can be written as the conjugation of the first inequality: 

K 1 + K,-L>«&~V-L)> V>- K M'- K * . (60, 

Now we write down the master rate equation for the probabilities p(L, z). We assume that it is a good approximation 
at large TV to replace sums by integrals, and write 

p(L, z) = J dK t J dK 2 G(K u K 2 \L) P {K u z)p(K 2 , z) - p(L, z) J C{L, K)p(K, z) (61) 

where the dot indicates differentiation with respect to z. The first term on the right-hand side of this equation gives 
the growth rate of L vortices from recombination processes as just described, and the second term reflects the loss of 
vortices by two-vortex interactions. The functions C, G are, in principle, to be determined from considerations of the 
full set of couplings in the scalar field theory of sec. IIIII 

Although we do not know these functions, they must obey certain conditions. There is one relation between them, 
following from conservation of total probability: 

0= J dLp(L,z)= J dKx [dK 2 p(K u z)p(K 2 ,z){ I ' dLG{K u K 2 \L) - C{K 1 ,K 2 )} (62) 

so that 

JdLG(K 1 ,K 2 \L) = C(K 1 ,K 2 ). (63) 

Furthermore, the essential conjugation symmetry p(L) = p(N — L) will be satisfied if the kernel G obeys 

G(K U K 2 \L) = G(N - K-l, N - K 2 \N - L) (64) 

with a similar equation for C . 

The probabilities we seek are those satisfying the equilibrium equation expressing independence of z (which we 
suppress) : 

p(L)=<yfdK 1 f dKtGiKuKilLMKtWK*), 7" 1 = f dKC(L,K)p(K). (65) 

At this point we drastically simplify the equations by taking both G and C to be constants Go, Cq within the regime 
of variables allowed by the conditions of Eqs. I|57|l . I|6U|) . Because these conditions are conjugation-symmetric, the 
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kernels continue to obey conjugation symmetry as in Eq. (|64|l . For suppressed processes not obeying these conditions 
we set G or G to zero. The decay function G is not subject to any particular restriction, so setting C to a constant 
means that the parameter 7 in the equilibrium equation 1)65(1 is independent of the p(K), by normalization of the 
probabilities. The constant Go must be 0(A _1 ), where A is the correlation length in z. The condition which preserves 
probabilities, Eq. (|63|l . then leads to Go = 0[(iVA) _1 ], because G is integrated over a range scaling like N, when its 
variables are in the scaling regime. This dependence of C, G on N means that the rate equation can be written solely 
in terms of the scaled probabilities p(x) defined in Eq. 1(38(1 . and we will do this below. 
The result is the approximate master equation 

P(L) = ^j Q dK 1 £ dK 2 {o[(N-L)- (N " Kl) ( N ~ K ^ ]+e(L-^^)}0[{L-K{)(L-K !l MK 1 )p{K 2 ). (66) 
The scaling form of this equation is 

p(x) = p(0) / dx\ I dx2\0[l - x ~ (\ — x\)(l — X2)] + 9{x — xiX2)\0[(x - x\)(x - X2)]p{xi)p(x2) (67) 



'0 Jo 

where we have made the replacement 



Go p(0) 

Go ~ N ' 



(68) 



That this is correct follows from Eq. H67fl by setting x = in the rate equation, which yields Eq. I|68|) . 

It is not our purpose in the present paper to study the rate equations in any great detail, largely because details 
of the solution will depend on details of the kernels G, G, which we do not know. However, there are certain features 
of the rate equations that we believe are robust: The behavior near zero, which is consistent with scaling, as we have 
already seen, and the behavior for flux L which is large compared to unity but still small compared to N/2. This is 
the regime of the modified DGA, already studied in sec. UVBl in which limits such as N or N/2 were replaced by 00, 
equivalent to replacing limits in the scaling variable x of ±1/2 by ±00. So in Eq. (|66|l . one sets the explicit As inside 
the integral to infinity (but not the N in the factor Go /Go), and in the scaling equation 1)68(1 one replaces the limit 
x = 1 by x = 00. After a little algebra, in the N — > 00 limit Eq. I(66|) becomes 

/•L />L />oo />oo 

p(L)= P (0) dKi p(K x )p(K 2 )+p(0) dK x 4Kip(Ki)p(Kz). (69) 



JO JL-Ki JL JL 

[In this equation there is a factor of A -1 hiding in p(0) and therefore a factor of N hiding in p(L).] 

For large L it is elementary to find the self-consistent behavior of the second term on the right-hand side of Eq. 1(69(1 , 
since in the integral only large values of Kx,K<2 contribute. The first term on the right has contributions from this 
regime but also from the regime where one or the other of K\ , K2 is small (if both are small, they do not contribute 
to large L). In this second regime we assume that p(K ) is constant and of order p(0) for K < p(0) _1 , and takes on 
its large- if behavior for K ~ L. After some analysis one finds that the self-consistent behavior at large L of p(L) is 
simply that of Eq. ((51|) and already given in Ref. |ll|, except that we can supply an overall factor: 

N 
p(L) = const— —2. (70) 

p{0)L 2 



or equivalently 



p(x) = const . (71) 



p(0) 



x 2 



where the constant factor is of order unity. One need do no analysis at all if one accepts that p{L) has a power-law 
behavior at large L, since L~ 2 is the only possible such behavior. 

This large-L behavior is quite acceptable in the scaling regime L ~ N, leading to densities uniformly of order 1/N. 
But there are certainly corrections for L in the non-scaling regime, since as already pointed out p(0) is finite. The 
simplest extrapolation that shows the asymptotic behavior at large L of Eq. 1(71) (I as well as the correct qualitative 
behavior at L = is the one already used phenomenologically in Eq. 1(47(1 : 

P{L) = *[{eN) 2 +L 2 Y P{X) = ,(e 2 +x 2 Y (72) 
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It follows that 

p(0) = - (73) 

and so the normalization of the large- L or large-x behavior in Eq. (|7U|I leads to the same parametric behavior in e as 
does the extrapolation formula given in Eq. J75J). Of course, we already know from sec. HVBl or if one accepts the 
DGA from (TJ, that this extrapolation leads to a string tension linear in k, at least for k <C N. 

It is simply not possible to argue for or against the extrapolation formula of Eq. I|72|) without knowing more than 
we know about the kernels C, G of the master equation. Let us, for example, modify the master equation (|67|l by 
replacing Go by a function G(L). In principle, the idea is that we are given G(L) and try to solve Eq. (|67|l for p(L). 
The solution is not obvious; it requires solving a non-linear integral equation. But one can think of the problem the 
other way around: Given p(L), to what kernel G(L) does it correspond? This question is answered by doing some 
integrals. If we use our guessed extrapolation in Eq. I|72|) we find that G(L) is a smooth and slowly-varying function 
approaching constant values as L — > 0, oo, no more or less plausible than assuming it is a strict constant. 

One may wonder how to recover conjugation symmetry from a non-conjugation-symmetric result such as the ex- 
trapolation of Eq. 1|72|) . One could make the substitution 

9 N 2 r ,2nL^ , s 

L ^2^ [1 " C0S( — )] - (74) 

The right-hand side above is indeed L 2 in the regime of interest to us: L <C N/2. Our extrapolation equation then 
becomes 

const 4ir 2 e 

PiL) = [ —\2ne) 2 + 2-2cos( 2 -^y (?5) 

To no one's surprise, this is precisely of the GO form given earlier in Eq. (|45|l if one identifies 



^(1)A 2 X , , 

27re = 2sinh( K ' ), (76) 

which in the small-e limit gives our previous result. 

VI. SUMMARY AND CONCLUSIONS 

Starting from the effective d = 3 gauge-theory action which yields center vortices and nexuses as solitons, we 
construct an TV-component scalar field theory. The components of this field theory are the unit-flux constituents 
(whose d = 2 cross-sections in a test plane we call atoms) of center vortices having higher flux J, whose d = 2 cross- 
sections we call molecules. Like constituents repel, and unlike constituents attract, with an attraction \/N times 
the repulsion. This leads to a condensate of merging, dissociating, and recombining vortices, and an approximate 
qualitative picture of vortex areal densities, based on studying these processes in a test plane that is pierced by 
vortices. This study leads in turn to a master rate equation, whose general structure we can foresee but whose 
detailed quantitative properties remain to be determined in future work. 

We expect that the basic dynamics of the d — 2 interaction of our atoms and molecules is the same in d — 4 and 
d = 3, and the master rate equation will have only minor quantitative differences in these two dimensions. 

We also give a rough and tentative first investigation of the rate equation, based on an approach to large- N dynamics 
which in some sense generalizes the usual implementation of the DGA. This generalization makes use of the work of 
GO, who propose a relation between fc-string tensions and areal densities (or equivalently probabilities) that is not at 
all restricted by the DGA. Its utility is that we can to a large extent bypass the complications of conjugation symmetry. 
At large N we replace discrete Fourier transforms by integrals with infinite limits, and show that the validity of this 
substitution defines a modified form of the DGA. These integrals preserve the scaling property p(L) = (\/N)p(L/N). 
Conjugation symmetry [p(L) — p(N — L), etc] is violated, but replaced by an equivalent condition that requires both 
p(L) and cr(k) to be even functions of their argument. We argue that this replacement of sums by infinite integrals 
violates no qualitative feature of the fc-string tension problem, provided that the scaling argument L/N is small 
compared to unity, and that k <C N/2. 

The rate equation, using some simplified kernels, has a solution p(L) that obeys scaling; which therefore has a 
probability of no vortex p(0) that is 0(1/N); and has an asymptotic behavior NL~ 2 (as proposed earlier in [ll| on 
phenomenlogical grounds). The simplest extrapolation of our asymptotic results to all fluxes is precisely the form 
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taken by the GO phenomenological form when L is large but obeys L <C N/2, and when a parameter e defining the 
modified DGA is small. This extrapolation is such that the fc-string tension is linear in k for k <C N. In this paper 
we leave entirely open the closely-related questions of the behavior which restores conjugation symmetry (requiring 
understanding the behavior for k > N/2), and of non-leading corrections to the large- N limit. 

Work is in progress to refine the drastic simplifications of the present paper, by calculating (from the gauge-theory 
effective action) those solitonic actions necessary to compute the scalar field theory coupling constants which describe 
merging, dissociating, and recombining vortices. Then one can hope to write down a more accurate master rate 
equation, properly including such effects as conjugation symmetry, and solve it numerically. 

Even without the elaborate work of calculating all the solitonic actions, it would be valuable to carry out computer 
simulations of the dynamics of atoms and molecules subject to the simple rules of attraction and repulsion used in this 
paper, incorporating all n-body processes of merger, dissociation, and recombination as well as conjugation symmetry 
and mod N flux conservation. 

It would also be of great value to study lattice simulations for moderately-large values of N, say, N = 4, 5, 6, and 
look for center vortices of various fluxes undergoing merging, dissociating, and recombining. In this way one could 
make a connection between direct simulations of A:-string tensions and the behavior of the vortex condensate. 

APPENDIX A: PROPERTIES OF VORTEX FLUX MATRICES 

We introduce the N x N traceless diagonal matrices 

Qj = diag( — , —,...—, -H , — , ...) i = l,2,...N (Al) 

N N N N N 

with the -1 in the ith position. Although N matrices are defined, their sum is zero, as is easily checked: 

A? 



Y^ Qi = 0. (A2) 



The independent Qi are just linear combinations of the weight vectors pi for the fundamental representation of SU(N) . 
The Qi obey 

Tr(Q t Q,) = Sa - ^ (A3) 

which is twice the scalar product jXi ■ p,j . 

Of course, any Qi is related to any other by a similarity transformation Qi — 5[ij]Qj5[ij] _1 . The S[ij] and their 
products form a representation of the permutation group, which means that various group orthogonality properties 
hold in summing over the vortex collective coordinates represented by the permutation group. 

The flux associated with each Qi is 2tt/N, in the sense that 

exp[2niQj] = exp[27ri/iV] (A4) 

for all j. Equation (|A2I) simply means that flux is conserved mod N. 
Next we define flux matrices for higher flux via 

F {J} = Q ll + ...Q lJ (A5) 

where no two indices are the same. If R of the indices occur at least twice each, we use the notation Fs j. m . Of course, 
to specify the detailed structure of the overlap we must write the index set in Eq. (|A5|) explicitly: {J} = {i\ . . . ij}. 
For example, the set {J} = {111223334567} has an overlap R — 3, but simply specifying R does not specify the exact 
overlap. In practice (see below) vortices with a given flux J and non-zero overlap always have higher action than for 
the same flux and no overlap, so in this paper we will restrict overlapped index sets to those in which an index in 
the overlap set {R} occurs exactly twice. The example above is not allowed, then, but {J} — {1122334567} is; it has 
overlap R — 3. Obviously, this restriction is easy to remove, but the consequent notation becomes cumbersome. 

The collective coordinates for -F{j} are permutations of the matrix indices — the same permutation for each Qi in 
the sum. In particular, if the set {</} is held fixed, and one sums over all permutations of the indices in another set 
{K}, there is orthogonality: 

J2 Tr[F {K} F {J} } = 0. (A6) 

permK 
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Note that there is a close connection between center vortices of flux {J} in the fundamental representation of 
SU(N) and elementary vortices (unit flux) in the J (totally antisymmetric) representation of this group. Every state 
of the J representation can be labeled by indices running from 1 to N as |ii . . . ij), with no two indices alike, and ket 
vectors with the same index set but in different order differ by the sign of the permutation required to go from one 
order to the other. For any diagonal matrix D — diag(Z?i, . . . I? at) in the fundamental representation, the action of D 
in the J representation on a ket is 

D\h...i J ) = (D il +...D ij )\i 1 ...ij). (A7) 

We have the equation 

ck TrjF ^ = ck TriF ^ ^ 

where Cj = J(N - J)(N + 1)/(2JV) is the quadratic Casimir of the J representation, and Dj = N\/[J\(N - J)!] is 
the dimension of this representation (note that J = 1 is the fundamental representation). On the left-hand side of 
Eq. (|A8(I traces are in the J representation; on the right-hand side, in the fundamental representation. 

So far the flux labels J are positive numbers. From Eq. IJA2JI it follows that if {J} and {iV — J} have no overlap 

F{j} + F{n-.j} = 0, -F{-,7} = ~F{j} = F{n-j}- (A9) 

One sees that 

exjp[2mF { ± J} ] = exp[±2mJN}. (A10) 

The master formula from which all other quadratic trace formulas follow is 

— IK 

Tr(F {J} F {K} ) = ^^ + R (All) 

if {J} and {K} have overlap (intersection) {R}- In this equation, both J and K are taken as positive numbers, so 
that if one uses this equation for a negative flux F{_,/}, one replaces this by — F{ J}. Or one comes to the same thing 
by using Fr N _n in place of — F{ J}, but then one must replace the overlap by R = K — R. 
The action of a center vortex {J} with no overlap is found by setting {J} — {K}, yielding 

^U = ■ M W A ^ (AID 

When there is overlap (as restricted above) we write the set {J;R} as the sum of two sets with no internal overlap, 
but of course the two sets overlap each other: {J; R} = { J — R} + {R}. One then finds the action for {J; R} to be 

TrFl J;R} = ^L^± + 2R. (A13) 

APPENDIX B: OVERLAP PROBABILITY 

Consider the temporary merger of vortices {J} and {K}, which have no internal overlaps but which taken together 
have overlap {R}. An exam ple is {J} = {12345}, {K} = {45678}, {R} = {45}. The composite vortex {J + K; R} 
has higher action [Eq. i|A13(l ] than a composite with the same total flux but with {R} — {0}. We wish to know, given 
J, K, the distribution of overlap number R. Evidently, given a set of K numbers taken from 1 to N, with no repeats 
allowed, the probability that one of an independent set of J numbers is in the set {K} is K/N. This suggests that 
(just as in tossing coins with a probability K/N of getting heads) the average overlap in J trials is 

<*> = ™- (Bl) 

The complete probability distribution of R, which wc call P(R\J, K), is found easily. If the overlap is indeed R, it 
occurs a number of times equal to the number of ways of picking J numbers subject to the constraint that R of them 
are in the set {R}, that is, the number of ways to choose J — R out of N — K. This is to be multiplied by the number 
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of ways an overlap of R numbers can be imbedded in {K}, which is the number of ways of picking R out of K, and 
divided by the total number of ways of picking J out of N. The result is 

P(mTK) - J\K\(N - J)l(N - K)l 

^(H\J,K) - mR ^ N _j_ K + R y( K _ R y {J _ R y ■ W 

Using Stirling's formula, with the assumption N 3> J, K, R 3> 1 , leads to a binomial distribution that is equivalent 
to a Poisson distribution: 

p(r\j,k) = —-H-—(^f(^-J^y-R _> LrlE)R e -jK/N . ( B3) 

v ' ' ; R\{J-R)\ K N> v N ' RV N ' K ' 

Note that the probability of no overlap (R = 0) is exp[— JK/N], which is exponentially small in the scaling regime, 
and that the expectation value of R is indeed JK/N. 
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